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I review the properties of the three-dimensional Gross-Neveu model formulated with 
non-zero chemical potential and temperature, focussing on results obtained by lattice 
Monte Carlo simulation. 

<N . i. INTRODUCTION 
(N 

^ ■ This talk concerns simulations of fermionic systems at non-zero density, which don't 

suffer from the problems associated with gauge theories. To paraphrase the sub-title of the 
workshop, they could be described as "not-quite-so complex systems with a real action" ; 
however, as I propose to show, the physics involved is still rich. The two theories I will 
discuss are referred to as "Gross-Neveu" models Q, and are described by the following 
p ■ Lagrangian densities: 
<D ! 

C A = fa($ +m )fa-jL(fafaY (1) 
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^ i C B = fa(@ + m)fa - {fafaf - {^Af (2) 

In each case the index % runs over Nf flavors of fermion, assumed to be four- component 
Dirac spinors. In the chiral limit m — > 0, symmetries akin to chiral symmetries can be 
identified: 

A Z 2 : j 5 ipi ; A f-> -4>a 5 (3) 

B U(l) : tpi exp(ia-f 5 )ipi ; fa fa exp(iatj 5 ) (4) 

A version of the model with a global SU(2) L <g> SU(2) fl symmetry can also be formulated. 

I will focus on the case where the models (1,2) are formulated in three spacetime 
dimensions, where the following features hold: 

• For sufficiently strong coupling g 2 the models exhibit dynamical chiral symmetry 
breaking at zero temperature and density. 

• The spectrum of excitations contains both baryons and mesons, ie. the elementary 
fermions, and composite fermion - anti-fermion states. 
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• The models have an interacting continuum limit. 

• For model B, whose chiral symmetry is continuous, the spectrum contains Goldstone 
modes. 

• When formulated on a lattice, the models have real Euclidean actions even for 
chemical potential /i ^ 0, and hence can be simulated by standard Monte Carlo 
techniques. 

For all these reasons, the models (1,2) are useful toys for understanding the behaviour 
of strongly-interacting matter at high density, as I shall discuss. 

2. MEAN FIELD ANALYSIS 

It is considerably easier to treat the Gross-Neveu model, both analytically and numer- 
ically, if the Lagrangians (1,2) are bosonised by introducing auxiliary fields: 

N f 

C' A = + m + a)^i + ji° 2 



N f 

£'b = M$ + m + a + in-f^i + -±(a 2 + 7i 2 ) (5) 

■^9 



Since the scalar field a and pseudoscalar it only appear quadratically without derivative 
terms, it is trivial to integrate them out and recover an identical generating function to 
that derived from the original Lagrangian. In the chiral limit m = the fermion may still 
have a dynamical mass S = (a). We can solve for the expectation value £ self-consistently 
via the Gap Equation, given by the diagram of Fig. la: 

S = 9 2 I tr -T^ ( 6 ) 
J P vp + L 

We find a non-trivial solution E ^0, which breaks the chiral symmetries (3,4), if 
J_ 1_ 2A 

9^9 9 J \ ' ) 

g 2 g l c vr 2 

where A is the UV cutoff on the momentum integral in (0). Note that £ — ► as g 2 /* g 2 . 
This result is exact up to corrections of 0(1/Nf), and it turns out that the self-consistent, 
or mean- field, approach coincides with the leading order of a 1/Nf expansion. 

To the same leading order in 1/Nf there is a correction to the scalar propagator (equal 
to g 2 /Nf at tree level) from the bubble diagram of Fig. lb. Remarkably, the linear 
divergence in this diagram is cancelled by the divergence in the definition of 1/g 2 , leading 
to a closed- form expression which is finite when expressed in terms of S: 

D ^ k ) = 7~7PY - (8) 

K> Nf(k 2 + 4Z2) tan- 1 (^) 

For the 7r field of model B, the factor (k 2 + 4S 2 ) is replaced by k 2 . We can examine the 
behaviour of (§) in two limits. In the IR regime k E, we find 

k 2 + \ 2E) 2 ; D M0C^, (9) 
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(a) 



(b) 



Figure 1. Leading order diagrams in the GN model 





Figure 2. 0(1/Nf) radiative corrections 



ie. the scalar resembles a weakly-bound fermion - anti-fermion composite whilst the 
pseudoscalar is a Goldstone mode. In the opposite UV limit k ^> £, we have 

1 



D a ,ir(k) oc 



(10) 



Thus the UV behaviour of the propagators is softer than would be expected for an aux- 
iliary field. Therefore higher order diagrams, such as those of Fig. |||, are less divergent 
than might be expected by naive power-counting: the conseqence is that the divergences 
can all be absorbed by rescaling the parameters of the original Lagrangian (||), and the 
1 / Nf expansion is exactly renormalisable . 

The transition between chirally symmetric and broken phases at g 2 = g 2 defines an 
ultra-violet fixed point of the renormalisation group. It is characterised by non-Gaussian 
values for the critical exponents: 



P 



1 = V 



5 = 2. 



Corrections to these values are 0(1/Nf) and calculable Indeed, they are currently 

known to 0{1/N 2 ) ||, and when extrapolated to small values of Nf are supported by 
Monte Carlo estimates ||. The continuum limit g 2 —>■ g 2 may be taken in either phase; 
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Figure 3. Mean field solution for E(/i, T) 



since our interest is in the restoration of the dynamically broken symmetry at non-zero 
density, we shall always work with a value of g 2 corresponding to the broken phase S ^ 0, 
(i/ji/j) 7^ 0. The rationale will be to refer all calculated quantities to a mass scale given by 
the physical fermion mass; to leading order in 1/Nf this scale is precisely E [§]. 

The mean-field gap equation (||) is readily generalised to non-zero temperature T and 
chemical potential fi: 

— = 4T V / — - - (12) 

g 2 n hJ (27r)2((2n-l)7rr-i M )2+^ + E 2 (/X)T) - ^ ) 

The equation can be solved @, eliminating the coupling g 2 which sets the cutoff scale in 
favour of the physical scale Eq, the value of E at \x = T = 0: 



E 



-T 



ln(l + e T ) + ln(l + e T 



(13) 



The phase diagram of the model in the (/z, T) plane consists of a chirally broken phase 
separated from a symmetric phase by the critical line E(/z, T) = 0, ie: 

T 

2— ln(l + e~r). (14) 



1 



JL 

So S 

For zero chemical potential we thus identify a critical temperature 
So 



T, 



21n2 ! 



whereas for zero temperature we identify the behaviour 
E(/x) = T, Q 9{ji c - /i) with fi c = E . 



(15) 



(16) 
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The phase transition is everywhere second order except at the isolated point T = 0: this 
can be traced to the slope of the surface S(/i, T) diverging in an essentially singular way 
as T — > 0, /j, — > \i c . Solutions for £(//) for various values of T/T c are shown in Fig. |. 



3. LATTICE SIMULATIONS FOR [l ^ 

The first Monte Carlo simulation of a four-fermi model with /j / was for the two 
dimensional Gross- Neveu model 0. The method is easily extended to three dimensions, 
and I will survey results from both models A M and B |]TD|. In all cases the staggered 
fermion formulation, involving single-spin component Grassmann fields x-, X is used; due 
to species doubling we have the result that in three dimensions N flavors of staggered 
fermion correspond to Nf = 2N continuum flavors [11 
the presence of a chemical potential is xM.\i with 



The free fermion kinetic term in 



5 ij 

J l xy 2 



u=l,2 



(17) 



and the sign factor r) v (x) = (— l) a ' oH ^"-i. 

The lattice version of the interacting theory starts from the bosonised form of the action 
@. For technical reasons it is preferable to formulate the auxiliary scalar fields a and 
7r on the dual lattice sites x []12|| . For model A, the full fermion matrix including the 
interaction with the scalars is thus 



Mg = M% + ~8 xy 8v £ cu 



X 



<x,x> 



where < x, x > denotes the set of dual sites neighbouring x. Note that M is real and 
diagonal in the flavor index i running from 1 to N. Noting the chiral symmetry analogous 
to (^), valid for bare mass m = 0: 



X{x) i-> e(x)x{x) ; x(x) i-> e(x)x(x) with e(x) = (-1) 



XQ+X1+X2 



(19) 



we see that the full global symmetry of the lattice action is O(N) <g> Z 2 . In the continuum 
limit, we expect this symmetry to enlarge to U(2iV) <S> Z%. Now, the hybrid Monte Carlo 



algorithm |L3[ simulates fermions by integrating over bosonic "pseudofermion" fields, using 
a real bilinear action <p(M* M) _1 0. Since M is real, the pseudofermion fields can be taken 
to be real, and the resulting expression for the path integral on integrating over cf) is 



Z= / Y[da(x)VdetM tr M cxp 



N 



2<7 S 



(20) 



For N even, the path integral measure due to the fermions is thus precisely detM as 
required. 

For model B the discussion is slightly more subtle, because now the matrix M is 
complex: 



Mg = M% + -5 xy 5v X (a(x)+ie(x)<K(x)). 

<x,x> 



(21) 
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Figure 4. £ vs. p for 1/g 2 = 0.7 

The weight in the functional integral after integrating over complex pseudofermion fields 
is thus det(.M + o + ieix) det(.M + a — ien), and is explicitly real (because the only 
complex terms occur on the diagonal of M). However, the simulation describes N flavors 
of fermion with positive axial charge and N flavors with negative axial charge, and the 
global symmetry including (f|) is U(iV) ® U(iV) ®U(1). In the continuum limit this should 
enlarge to U(2iV) ® U(2iV) ® U(l), but it is important to note that the U(JV» ® U(l) 
symmetry of (2) is not recovered. Even for \x = 0, the action (2) leads to a complex 
measure, and cannot be simulated by this algorithm. 

Model A has been simulated using the partition function ( p0|) with Nf = 12 0. Since 
the chiral symmetry group is discrete, and since the diagonal elements of M are in general 
non-zero, it is possible to work directly in the chiral limit m — 0. Fig. f| shows £(//) 
evaluated on a 20 3 lattice (approximating T = 0) with coupling 1/g 2 = 0.7, corresponding 
to the system being in the broken phase at zero density with So = 0.43 (the bulk critical 
coupling 1/g 2 is estimated to be 0.98(3)). We see that the order parameter S is very 
weakly dependent on /i up to a critical value /z c ~ 0.394 ~ E , whereupon it drops to zero 
dramatically, suggesting a strong first order transition in accordance with the mean-field 
prediction. Also plotted is the fermion number density p, defined by 

p( /i ) = _I^£ = J_/y Um- 1 , + e^M" 1 A (22) 
V dp 2V\ ^ x,x+o x,x-oj /' v * 

which shows a sharp rise across the symmetry-restoring transition. It is tempting to 
interpret the p ^ signal for p < p c as evidence for a "nucleon liquid" , but this would be 
premature, since it is quite possibly a finite volume (ie. non-zero temperature) artifact. 
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Figure 5. X vs. pon asymmetric lattices 



Indeed, test simulations on a 12 3 system show that the transition is softened considerably 
on a smaller volume, with p increased for \i < fi c . It is also possible to study the approach 
to the continuum limit. Since fi couples to a conserved current, it should not undergo 
renormalisation, and hence \x c should scale like a physical mass, leading to the relation 

^ c0C \ ~^~ ' ( 23 ) 
Wc r) 

with v the correlation length critical exponent. From runs at 1/g 2 = 0.7, 0.75, 0.8 we find 
relation (^) verified with v = 1.05(10), in agreement with the leading order prediction 

We have also studied the transition at non-zero temperature by simulations on asym- 
metric lattices N T x 36 2 with 1/g 2 = 0.75. The results are shown in Fig. [5| The data show 
that the transition is softened, with \i c a decreasing function of temperature, in qualitative 
agreement with the mean field prediction. Estimates of the location of the phase transi- 
tion are plotted in Fig || together with the mean field prediction fll4|) , showing reasonable 
agreement; the discrepancies can probably be accounted for by 0{1/Nf) corrections. 

There are two issues raised by the lattice simulations that are worth discussing. First, 
what is the order of the transition for T > 0? The isolated nature of the first order point 
predicted by mean field theory seems unlikely to persist once quantum effects are included 
in the calculation. Indeed, the simulations showing a two-state signal characteristic of the 
first order transition were performed on lattices of finite extent, and by implication at 
a small but non-zero temperature. The most likely scenario is for the line of first order 
transitions to extend into the (/x, T) plane, to end at a tricritical point separating it from 
a line of second order transitions continuing to the point (0, T c ). This behaviour has been 
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Figure 6. Phase diagram in the (/x, T) plane 



verified by calculation of 0(1/Nf) corrections in the two dimensional model [fn}| , although 
ironically, the thermodynamics in this case are not described by the 1/Nj expansion due 
to condensation of topological defects in the resulting one dimensional system fT5|j M. 
Numerical studies with higher precision, or explicit calculations of 0(1/Nf) corrections 
in three dimensions would be useful to clarify the situation further. 

Secondly, what is the nature of the pure thermal phase transition, ie. along the fi = 
axis? The mean field solution (|13|) predicts a transition with Gaussian exponents (3 = v = 
2, 8 — 3, 7 = 1, rj = 0, and indeed, early simulations 0] on a 12 x 36 2 system with Nf = 12 
and with m^0, estimated 5 = 3.0(1) in apparent confirmation. This led to speculation 
[[HJ that models in which a composite field (in this case a bound fermion - anti-fermion 
pair) condenses violate the conventional wisdom |17 |, which predicts that the thermal 
transition has the exponents of the dimensionally reduced ferromagnetic model with the 
same global symmetry, in this case the two dimensional Ising model. More recently there 
has been an analytic study [|TB[ using both dimensional reduction and subsequent linked 
cluster expansion techniques, and a more accurate numerical study Ul9| : both favour the 
2d Ising model exponents for Nf < oo. 

Next I turn to numerical studies of Model B [|10j, which has a continuous U(l) chiral 
symmetry and hence a Goldstone pion in the broken phase. Here all simulations were 
performed on a 16 3 system with bare mass m = 0.01, but this time with differing numbers 
of flavors, ie. Nf = 4 and 12 in continuum normalisation. Just as for the Z2 model, it was 
found that a sharp chiral symmetry restoring transition occurs at a critical /i c ~ O(So), 
and that bulk expectation values of number density p and energy density were in good 
qualitative agreement with mean field theory, the greater discrepancy occuring as expected 
for Nf = 4. More interesting results come from studies of the spectrum; the fermion 
mass ruf is measured by studying decay of the timeslice propagtor G + (t) in the timelike 
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Figure 7. Spectrum results for the U(l) model 



direction (a ± sign denotes that G is calculated with a factor e ±At on forward-pointing 
links), while the decay of the auxiliary field correlator yields the pion mass m^. Results 
are shown in Fig. [?| for Nf = 12, 1/g 2 = 0.5 (ie. small quantum corrections, large lattice 
spacing), and Nf = 4, 1/g 2 = 0.7 (larger quantum effects, closer to the continuum limit). 

What is observed is a clear separation between the scale defined by the critical chemical 
potential /i c ~ m/(/i = 0), and the scale m n /2. This is significant, since it implies that the 
Gross- Neveu model does not suffer from the early onset of symmetry restoration observed 
in simulations of both quenched ||2(J and full |21j QCD. Instead we find rrif(fi) ~ m/(0)— //, 



and m n roughly constant for /i < /i c , and rising sharply for /i > fx c , in accordance with 
naive expectation. Why should this be? The answer seems to be that the Goldstone 
mechanism in these models is realised through a pole in the pseudoscalar channel formed 
from disconnected diagrams (ie. sum of bubbles), in accord with the 1/Nf expansion. 
Indeed, the connected (ie. flavor non-singlet) diagram, formed from the expectation 
(G + (t)G_(t)), appears numerically much less important, yielding a bound state mass of 
roughly 2mj(0) p2[ . If the connected diagram had been numerically dominant, as is the 
case in gauge theories, then the same diagram would have predicted a light state formed 
between two fermions of opposite axial charge (recall the discussion below eqn. (|2TD); 
since this state carries net fermion number, it would be equivalent to a baryonic pion, 
and hence promote a premature restoration of the symmetry as fi increases |23| . 



4. SUMMARY AND OUTLOOK 



The Gross-Neveu model presents an "easy" problem, in the sense that it has a well- 
defined continuum limit, and can be explored in a controlled fashion without introducing 
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ad hoc interactions or form factors. Nonetheless, the physics it contains for /x 7^ is rich, 
and is possibly not too dissimilar, at least qualitatively, from that of QCD. Issues such 
as the nature of the thermal transition, the existence of a tri critical point, and also a 
possible nucleon liquid-vapour phase transition (expected for /x < /x c for sufficiently low 
temperatures) can all be addressed by more thorough simulations than those reviewed 
here. Moreover, quantitative information is also available in principle from the 1/Nf 
expansion, known to be accurate for Nf as small as 2 ||. In addition, the model might 
also serve as a useful warm-up exercise for other analytic methods such as the Schwinger- 
Dyson or exact renormalisation group approaches. Despite its apparent easiness, however, 
the model also presents a non-trivial test for numerical algorithms, such as the Glasgow 



reweighting method ETJ, designed to treat systems with complex integration measures. 



One insight to come out of the study p2| is that the Gross-Neveu model appears to 



be simulable by conventional means because the Goldstone mechanism is realised in a 
fundamentally different way to the case in gauge theories, reflecting in turn the fact that 
the global symmetries of the two problems are different for general Nf. 

Let me finish by advertising the fact that there is also a "difficult" problem in three 
dimensional fermion physics, namely the Thirring model, with the following Lagrangian 
density: 
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C T hirring = 1pi($ + m)4>i + — — (^iT^i) • (24) 



2N f 



This model has been studied by 1/Nf expansion ||24j| , Schwinger- Dyson equations J25 
and numerical simulation | 26f . I believe it to be the simplest fermionic system for which 



a numerical solution is necessary. Its essential features are: 

• It is non-perturbative in 1/Nf. 

• It exhibits dynamical chiral symmetry breaking and a UV fixed point of the renor- 
malisation group for Nf < 4. 

• The chiral symmetry group is U(iV/), similar to QCD. 

• For /x 7^ the action is complex. 

In other words, the model exhibits all the essential physics except confinement. Perhaps 
it would provide a suitable intermediate challenge while we wait for the means to solve 
the real problem. 
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